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ABSTRACT 

Estimation of biochemical constituents in vegetation such as 
lignin, cellulose, starch, sugar and protein by remote sensing methods is 
an important goal in ecological research. The spectral reflectances of 
dried leaves exhibit diagnostic absorption features which can be used 
to estimate the abundance of important constituents. Lignin and 
nitrogen concentrations have been obtained from canopies by use of 
imaging spectrometry and multiple linear regression techniques. The 
difficulty in identifying individual spectra of leaf constituents in the 
region beyond 1 pm is that liquid water contained in the leaf dominates 
the spectral reflectance of leaves in this region. By use of spectrum 
matching techniques, originally used to quantify whole column water 
abundance in the atmosphere and equivalent liquid water thickness in 
leaves, we have been able to remove the liquid water contribution to 
the spectrum. The residual spectra resemble spectra for cellulose in the 
1.1 pm region, lignin in the 1.7 pm region, and starch in the 2.0-23 pm 
region. In the entire 1.0-23 pm region each of the major constituents 
contributes to the spectrum. Quantitative estimates will require using 
unmixing techniques on the residual spectra. 


INTRODUCTION 

The determination of leaf chemistry from spectral remote sensing of vegetation 
canopies is a major goal in terrestrial ecology because it would make possible detailed study 
of biogeochemical processes such as productivity, decomposition and nutrient turnover rates 
on regional and global scales not otherwise accessible. So far, spectral reflectance 
measurements in the 0.4-25 pm region have been applied to dry ground leaf matter and the 
results have yielded good correlation between spectral reflectance in individual narrow (10 
nm) bands and leaf components such as nitrogen and lignin [1,2], In fresh, green leaves the 
reflectance spectra beyond 1.0 pm are dominated by ligand water absorption both within 
and without the strong 1.4 and 1.9 pm absorption features. The dominant water absorption 
to a great extend obscures the diagnostic spectral features of the leaf biochemical 
components. In this paper we describe a technique to subtract the water contribution by 
modeling and obtain residual spectra that exhibit characteristic leaf-component spectra. 

BACKGROUND 

Changes in ecosystem processes such as productivity and decomposition may be 
expressed in the canopy foliage chemistry as a result of altered carbon allocation patterns, 
metabolic processes and nutrient availability. Partitioning of carbon and nutrient resources 
between shoots and roots is reflected in the relative concentrations of carbohydrates and 
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nitrogenous compounds in leaf tissue [3,4]. Environmental stress, including change in nutrient 
availability, may alter resource partitioning within the plant. Increases in nitrogen 
supply frequently results in increased foliar nitrogen concentration [5,6,7], whereas 
reduction in availability promotes increased secondary wall thickening and lignification 
[8,9]. Storage carbohydrates such as starch, and defensive compounds, such as polyphenols 
and fibers, vary predictably with resource availability and heibivory pressures [10,11,12]. 

The quantity and organic chemical composition of litter produced by the canopy are 
important controlling factors in the processes of decomposition, mineralization and 
nitrification [13,14,15,16]. The recalcitrant nature of lignin modulates the rate at which 
litter decomposes and, hence, the release of nutrients bound in the dead tissue. Significant 
changes in foliar nitrogenrlignin ratios may indicate corresponding changes in 
decomposition rates affecting nutrient cycling and trace gas fluxes [17,18,19]. 

The assessment of canopy biochemistry using remote sensing would provide 
knowledge of the spatial extent and variation of carbon/ nutrient sources and sinks crucial to 
understanding gas exchange between vegetation and the atmosphere. Research in 
analytical chemistry has demonstrated that concentrations of constituents within organic 
mixtures can be evaluated from near infrared reflectance spectra of those mixtures [20,21]. 
Recent work in imaging spectrometry indicates that high spectral resolution data contain 
more information on subtle spectral features than do broad-band data. Canopy lignin 
concentrations were successfully estimated using NASA’s Airborne Imaging Spectrometer 
(AIS) and were subsequently used to derive images of annual nitrogen mineralization rates 
in a series of Wisconsin forest ecosystems [22,23]. 

In the AIS study, stepwise linear regression techniques were used to determine the 
single spectral band or bands that best correlated with measured lignin content [22,23]. The 
technique does not take advantage of the spectral shape of the lignin reflectance feature as 
a diagnostic tool nor the increased discrimination among other leaf components with 
overlapping spectral features 

LABORATORY MEASUREMENTS 

High spectral reflectance measurements of stacks of five fresh oak leaves were 
made in a Beckman 5270 spectrophoto meter utilizing an integrating sphere attachment. 
Reflectance against a Halon standard was acquired at 5 nm intervals between 0.8 and 2.5 
|im. The specimens were weighed before and after each 45 minute run to assess dessication. 
The weight loss was less than 10%. Spectra for fresh and dried leaves were measured.The 
water model spectrum was obtained by wetting a sample holder containing 100-250 pm 
diameter glass beads that act as inert scattering centers because of refraction between glass 
and water. The beads mimic the effects of leaf cell structure that contains numerous 
refractive index boundaries. Spectra of dry glass beads show no spectral absorption features 
and, therefore, the wetted-bead spectrum yields the pure water component of the leaf 
spectrum. 

SPECTRAL CURVE FITTING 


Figure 1 shows measured reflectance spectra of fresh oak leaves (/^(A)), dry oak 
leaves and glass beads mixed with liquid water Because the glass 

beads alone have no absorption bands in the 1. 0-2.5 |im region, the R WM (X ) contains only 

liquid water absorption features centered at approximately 1.18, 1.43, and 1.92 |im. In order 
to model the observed spectra properly, the absorption and scattering properties of the 
sample must be known and rigorous multiple scattering calculations are required. This kind 
of modeling is difficult in practice. We have developed a spectral curve fitting technique to 
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model reflectance spectra of fresh leaves. The procedures are as follows: 

1. The approximate absorption coefficients of dry leaves (AT 0 (A)) and liquid water 

{K W (X )) are derived from the reflectance spectra R DM {X ) and R WU (X ) according to the 

Eqs. K D (X ) = — log(/? OAf ( A )) and K W {X ) = — log(/?^ w (A )). These reflectance-based 

absorption coefficients have already included the multiple scattering effects. They are not 
the same as the commonly used absorption coefficients, which are based on transmittance 
measurements. 

2. The fresh leaf spectrum )) is calculated according to 

R FC {X ) = (A + fiA)x exp {-\C W K W {X ) + C D K D (X )]} (l) 

where A, B, C w and C D are adjustable parameters. The term (/t + BX ) represents the 
background level of the calculated spectrum. 

3. The parameters in Eq. (1) are adjusted so that the sum of the squared differences between 

the observed spectrum R FM and the calculated spectrum R FC in a given spectral region is 
minimized. The processes of parameter adjustment and minimization are automated by 
incorporating our spectral calculation program into a non-linear least squares fitting 
program [24]. 

For purposes of discussion, we show the fitting procedure only for the 15-1.7 pm region 
where lignin is an important constituent. Similar results have been obtained in the 1.1 pm 
region where cellulose is influential and in the 2.0-2.3 pm region where the starch 
component is important. 

Fig. 2 shows an example of spectral curve fitting. Both K w and K D are included in the 
spectral calculations, and an excellent fit between the calculated and the observed spectra 

is obtained. Fig. 3 shows a similar fit, but excludes K D in the spectral calculation (or 

equivalently, fixes C D to 0), and a poor fit is obtained. Fig. 3 also shows the residual 
spectrum, which resembles the dry leaf spectrum. Therefore, the reflectance spectra of 
fresh leaves in the 1.6 pm region are controlled by absorptions of both the liquid water in 
the leaves and the dry leaf components. This contradicts conclusions of Allen et al. [25] that 
a sheet of liquid water can completely account for the absorption spectrum of a single leaf in 
the 1.4- to 25-pm spectral range. 



Fig. 1. Reflectance spectra of fresh and dried oak leaves and glass beads plus water. 
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Fig. 2. Fit of the fresh leaf spectrum with the dry leaf and water spectral reflectance 
components. 



Fig. 3. Fresh leaf spectrum fitted with water spectral reflectance only. Below is the 
difference between the two spectra. 

RESULTS AND CONCLUSIONS 

Reflectance spectra of fresh green leaves are composed of dry leaf spectral 
components and water spectra. Therefore, by the proper application of model water 
spectra, it is possible to derive a residual spectrum that resembles the dry leaf spectrum. 
Further decomposition of the residual spectra into individual components should be 
possible using unmixing techniques [26]. 

The curve fitting technique will be applied to canopy spectra obtained from the 
Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) data, and it is expected that 
the equivalent water thickness derived from the 0.98 pm liquid water feature [27] can be 
used in the modeling. 
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